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On  The  Calculating  Models  of  Permanent  Magnets 

Sun  Yushi ® 

© 


© 


(Nanjing  Aeronautical  Institute) 

[Abstract]  Two  improvements  were  proposed  for  two  calculating 
models  for  permanent  magnets  (the  scalar  potential  model  and 
the  vector  potential  model)  in  this  paper.  Volumetric  density 
was  replaced  by  an  appropriate  hypothetical  surface  density 
of  magnetic  monopole  (or  surface  bound  current  density)  which 
did  not  vary  with  the  operating  point.  The  definition  of  mag¬ 
netic  reluctivity  was  correspondingly  modified  in  the  vector 
potential  model  to  simplify  the  calculation  and  computer  pro¬ 
gram.  The  improved  models  have  been  proven  through  computa¬ 
tion  and  experiments. 

I .  Introduction 


The  widely  used  calculating  models  for  permanent  magnets 
in  engineering  include  the  scalar  potential  model ^  (Model  I 

[1'21  (Model  II 


hereafter)  and  the  vector  potential  model 
hereafter) . 

According  to  Model  I,  a  permanent  magnet  is  considered  to  be 
comprised  of  distributed  hypothetic  magnetic  monopoles  (the 
volumetric  magnetic  monopole  density  is  and  various 

anistropic  magnetic  conductors.  Therefore,  a  permanent  mag¬ 
netic  field  is  expressed  by  the  following  quasi-Poisson 
equation . 


divugradu  *  -P 


mO 


(1) 


and 


Pm0  *  -“odivMo 


(2) 


where  U  is  the  scalar  magnetic  potential,  Hg  is  the  magnetic 
permitivity  in  vacuum,  and  Mg  is  the  magnetization  vector  in 
the  permanent  magnetic  when  the  magnetic  field  strength  H  is 
© 


zero. 
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3  For  ease  of  expression,  rectangular  coordinates  were  used 
throughout  this  work.  It  was  also  assumed  that  the  original 
permament  magnet  and  magnetizing  directions  were  along  the 
Z-axis. 


The  definitions  of  BQ^  and  are  shown  in  Figure  1.  *  is 


the  magnetic  permitivity  matrix  of  the  permanent  magnet: 


Figure  1 . 


Geometric  Significance  of  Several  Magnetic  Conduc 
tivity  (or  Magnetic  Reluctivity) 
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(B) 

Furthermore 

(3) 

(4) 

According  to  Model  II,  a  permanent  magnet  is  considered  to 
be  comprised  of  bound  currents  (the  volumetric  current  density 
is  JmQ)  and  various  anisotropic  magnetic  conductors.  There¬ 
fore,  a  permanent  magnetic  field  (when  a  macroscopic  current 
density  J  is  present)  satisfies  the  following  equation 

rot<vr9tA)-Jme+j 

(5) 

Jmt—  0,rotO'Mi> 

(6) 

where  A  is  the  vector  magnetic  potential  and  v  is  the 
magnetic  reluctivity  matrix,  and 

apparent 

!*  >,  0  O' 

*m-  o  v,  o  |,  — L_, 

-  C  0  V.  J 

**  9 

(C) 

In  the  two  models  mentioned  above,  u.  or v  „  cannot  be 

^3  q 

found  from  experimental  data.  Furthermore,  their  selection 
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is  not  consistent.  Frequently,  a  certain  parameter  in  the 
magnetization  direction  of  the  permanent  magnet  is  made  to 
be  Uq. 

The  author  made  the  following  modifications  through 
supplementing  derivation  and  reorganization  of  equations  (1) 
and  (2) .  The  expressions  for  surface  monopole  density  and 
bound  surface  current  were  supplemented.  In  Model  II,  the 
bound  volumetric  current  could  be  avoided  in  ordinary  conditions 
through  an  appropriate  modification  of  magnetic  reluctivity. 
Furthermore,  bound  surface  current  was  made  to  be  independent  of 
the  non-linear  coefficient  matrix  to  simplify  the  calculating 
process.  In  addition,  objections  were  raised  against  the 
determination  of  magnetic  permitivity  u  (or  magnetic  reluc- 
tivity)  perpendicular  to  the  direction  of  magnetization  in  the 
original  model.  It  was  pointed  out  that  must  be  obtained 
experimentally . 

II.  Modification  of  Model  II 

The  characteristic  equations  (3)  and  (4)  for  Figure  1 
(a,b,c)  can  be  expressed  by  the  following  matrix  equation: 

(7) 

By  taking  the  curl  on  both  sides  of  eq.  (7) ,  we  got  equations 
(5)  and  (6).  On  this  basis,  the  boundary  conditions  of  the 
medium  with  macroscopic  surface  current  were  compared  to  those 
without  such  current (surface  bound  current  density  is  considered 
as  macroscopic  surface  current  density  in  Model  II) .  An  ex¬ 
pression  for  the  surface  bound  current  density  jmQ  could  be 
derived  from  the  tangential  component  of  H. 


X  n  ^  0 

where  n  is  a  unit  vector  in  the  normal  direction  along  the 

surface  of  the  permanent  magnet.  From  equations  (6)  and  (8) 

one  can  see  that  J  „  and  are  both  directly  related  to  the 

mo  wmo 

vector  ug(uMo). 
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If  uQ  (vMQ)  a  constant/  then  JmQ  is  always  zero  ^  and  jmQ 
is  a  constant.  However,  »  usually  varies  with  the  operating 
point  for  a  permanent  magnet.  Therefore,  uQ  (vMQ)  cannot  be 
a  constant.  Hence,  in  Model  II  which  is  expressed  by  equa¬ 
tions  (5) ,  (6)  and  (8) .  It  is  inavoidable  to  have  JmQ  present 
and  jm0  varying  with  the  operating  point. 

There  were  enough  reasons  to  express  the  characteristics 
of  a  permanent  magnet  [see  Figure  1  (d,  e)]  by  another  matrix 
equation  i . e . , 

where  is  the  equivalent  magnetic  reluctivity  matrix  of  the 
permanent  magnet.  Furthermore, 


r  »•».  o  o  • 

*•»*!  0  Vm,  0 

L  0  0  vx,w 

“  ( H  i  —  H  i  j  y/B  i 

Hq*  is  the  equivalent  goercivity  vector  and 


r  i 

m :  ~~Ht,  !  (D) 

L  -H„  J 

Obviously,  equation  (9)  is  equivalent  to  (7) .  Furthermore, 
under  a  linear  condition,  » v  and  Hq  *  uq(vMq).  After 
treating  equation  (7)  with  the  same  treatment  as  for  equation 

4  This  conclusion  is  applicable  to  a  permanent  magnet  magnet¬ 
ized  in  the  z  direction  and  along  a  radial  direction. 
However,  it  is  not  suited  for  circumferential  magnetization. 
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(9),  the  overall  expressions  for  the  improved  Model  II  were: 


A)=Jmh  +  J  1 
■lmn  =  :c:H 

1  mu  —  Ht  X  n  J 


(11) 


From  the  magnetization  direction  it  is  possible  to  determine 
that  J  is  zero  in  most  cases  (such  as  magnetizing  along  z- 
axis  or  radially) .  It  is  a  constant  vector  only  when  it  is 
magnetized  circumferentially.  However,  j  ,j  is  invariantly  a 
constant  vector. 

From  the  equivalence  of  mathematical  expressions,  the 


selection  of  H^  and  is  arbitrary,  in  principle.  However, 


whether  the  choice  is  appropriate  will  affect  the  converging 
rate  of  the  computation.  For  example,  in  example  2  in  this 
paper,  when  the  permanent  magnet  is  samarium-praseodymium- 
cobalt  [its  characteristic  is  shown  in  Figure  1  (f ) ] ,  the  op¬ 
erating  region  of  the  permanent  magnetic  is  mainly  in  the 
linear  portion  of  the  curve.  Therefore,  the  fluctuation  of 


vHp  can  be  minimized  in  iterations  by  choosing  HQp  at  the 


intercept  of  the  extension  of  the  linear  section  of  the  char¬ 
acteristic  curve  with  the  horizontal  coordinate  obtaining  the 
fastest  convergence. 

III.  Supplement  to  Model  I 

From  equation  (7) ,  the  boundary  conditions  of  the  media 
with  a  hypothetic  surface  magnetic  monopole  density  were  com¬ 
pared  to  those  without  it.  The  expression  for  the  hypothetical 
surface  magnetic  monopole  density  on  a  permanent  magnet  could 
be  derived  by  using  the  normal  component  of  B: 


•/»  (12) 

From  equation  (2)  and  (12)  one  can  see  that  both  amQ  and  pmQ 
are  unrelated  to  n.  Therefore,  when  MQ  is  a  constant,  pmQ  in 
most  cases  is  invariantly  zero.  It  is  a  constant  only  when 
the  magnetization  is  along  a  radial  direction  Z.om0,  however, 
is  always  a  constant.  This  shows  that  in  most  cases  only  the 


effect  of  surface  magnetic  monopole  must  be  considered.  Volu¬ 
metric  magnetic  monopole  does  not  have  to  be  taken  into  account 
IV.  Magnetic  Permitivity  in  Perpendicular  Direction 

The  author  conducted  experiments  on  two  batches  of  speci¬ 
mens  made  primarily  of  AlNiCo-5  and  the  typical  data  and 
curves  are  shown  in  Table  1  and  Figure  2.  The  experimental 
results  explained  the  following  three  situations,  is  an 
independent  parameter  different  from  parameters  such  as  u  , 

icC 

Kr,  and  4^  which  are  determined  by  the  magnetization  direction. 

Within  the  range  of  H  <  0 . 5_H_ ,  y  _  is  nearly  a  constant .  When 

q  B  t  q 

is  not  large,  the  hysteresis  of  the  Bq"Hq  curve  is  very 
small  and  Mg^  can  approximately  be  considered  as  zero  for  a 
permanent  magnet.  is  not  apparently  related  to  the  magnetic 

state  in  the  magnetization  condition.  Therefore,  u  should  be 

4 

obtained  experimentally. 

Table  1 .  Average  Values  of  Important  Parameters  of  the  First 
Batch  Samples  (Seven)  in  Magnetization  and  Perpendi¬ 
cular  Directions 

3*C  i  fi, 

<A/a>  ( lO'*H/m) i  f  lu*  H/m).  ( I0"*H.  m) 


3r 

■  T ) 


Table  2.  Typical  B  -H  Characteristic  Curve  of  a  AlNiCo-5 

Specimen  (reproduced  according  to  the  curve  measured) 


■3, 


1.  scale  Mg  =  Q.15T/div.  M^  =  4.07  x  10^  A/m/dov. 

V.  Experimental  Proof 

Example  1.  A  square  cross-section  permanent  magnetic  bar 
was  used  to  prove  the  expression  for  surface  magnetic  monopole 
in  Model  I.  Furthermore,  the  calculated  curve  of  and  the 
measured  points  are  shown  in  Figure  3. 

The  finite  difference  method  was  used  in  theoretical 
calculation.  A  coarse  computation  was  carried  out  in  a 
relatively  large  area  (nodal  points  27  x  27  x  26  *  18954)  . 
Then,  a  fine  calculation  was  performed  in  a  correspondingly 
reduced  area  (nodal  points  26  x  26  x  25  *  16900) .  The 
demagnetization  and  u  used  in  the  computation  were  all 


actually  measured  numerical  values. 


Figure  3.  Calculated  Curves  and  Distribution  of  Measured 
Points  (B^  in  Example  1. 

The  measuring  instrument  was  a  Hall  effect  gaussmeter.  From 
Figure  3  one  can  see  that  the  theory  coincides  with  practice. 
Example2.  A  permanent  magnetic  ring  (made  of  YX-30)  in  the 
radial  direction  was  used  to  verify  the  improved  Model  II  (see 
Figure  4) .  For  ease  of  comparison,  calculated  results  obtained 
based  on  Model  I  are  also  given  in  Figure  4  (the  finite  differ¬ 
ence  method  was  used  for  Model  I  with  95  x  34  =  3230  nodal 
points) . 

Both  the  finite  difference  method  and  the  finite  element 
method  were  used  in  the  verification  process  (the  grid  division 
of  the  finite  difference  method  was  the  same  as  that  in  Model 
I,  and  the  nodal  points  were  164  for  the  finite  element  method) . 
The  results  are  basically  in  agreement  with  those  in  Model  I. 
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The  selection  of  Hrt  and  v„  is  shown  in 

Op  Hp 


In  addition,  the  permanent  magnetic  ring  was  used  to 
verify  the  accelerated  convergence  effect  of  the  improved 
Model  II.  The  permanent  magnet  shows  a  non-linear  character 
at  below  0.7T. 

Figure  1  (f) .  The  original  and  improved  models  were  used  to 

calculate  the  potential,  and  the  under  relaxation  iteration 

method  was  used  to  correct  the  non-linear  coefficients  v_  and 

P 

v  .  The  converging  situations  using  various  under-relaxation 
Hp  CO 

factors^are  listed  in  Table  2.  On  can  see  that  the  converging 
rate  is  greatly  improved  after  Model  II  was  improved. 
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Figure  4.  Comparison  of  Results  Obtained  Using  Finite  Differ¬ 
ence  Method  with  Model  II  to  Model  I  (BR  portion) 

1.  improved  Model  II  2.  Model  I 
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Table  2.  Comparison  of  Convergence  Under  Various  W 
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